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Properties of collective flow and pion production in intermediate-energy heavy-ion collisions with a 
relativistic quantum molecular dynamics model* 


Si-Na Wei! and Zhao-Qing Feng! t 
'School of Physics and Optoelectronics, South China University of Technology, Guangzhou 510640, China 


The relativistic mean-field approach was implemented in the Lanzhou quantum molecular dynamics trans- 
port model (LQMD.RMEF). Using the LQMD.RMF, the properties of collective flow and pion production were 
investigated systematically for nuclear reactions with various isospin asymmetries. The directed and elliptic 
flows of the LQMD.RMEF are able to describe the experimental data of STAR Collaboration. The directed flow 
difference between free neutrons and protons was associated with the stiffness of the symmetry energy, that is, 
a softer symmetry energy led to a larger flow difference. For various collision energies, the ratio between the 
a and m” yields increased with a decrease in the slope parameter of the symmetry energy. When the collision 
energy was 270 MeV/nucleon, the single ratio of the pion transverse momentum spectra also increased with de- 
creasing slope parameter of the symmetry energy in both nearly symmetric and neutron-rich systems. However, 
it remained difficult to determine the dependence of the double ratio on the stiffness of the symmetry energy. 
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I. INTRODUCTION 


The equation of state (EOS) of nuclear matter, which orig- 
inates from the interactions between nuclear matter, plays an 
important role in nuclear physics and astrophysics. Heavy- 
ion collisions, the properties of nuclei, and neutron stars 
(NSs) have been widely studied to extract the nuclear EOS. 
Because nuclear many-body problems are highly nonlinear 
and the EOS is not a directly observable quantity in experi- 
ments, there are still some uncertainties in the EOS despite 
great efforts [1-6]. For instance, the EOS extracted from 
the GW170817 event has uncertainties at high nuclear den- 
sities [2]. Although the EOS can be extracted from the prop- 
erties of NSs, the internal composition of NSs is still poorly 
understood. The core of an NS may contain exotic mate- 
rials such as hyperons, kaons, pions, and deconfined quark 
matter [7-12]. Heavy-ion collisions in terrestrial laboratories 
provide a unique opportunity to study both the EOS and ex- 
otic materials. 

Collective flows of heavy-ion collisions were proposed in 
the 1970s and first detected in experiments at Bevalac [13- 
16]. Because collective flows are associated with nucleon- 
nucleon interactions, nucleon-nucleon scattering, etc., collec- 
tive flows have been used to extract the nuclear EOS [1]. Col- 
lective flows are also helpful for understanding the phase tran- 
sition between hadronic and quark matter. Generally, when a 
phase transition between hadronic and quark matter occurs, 
the collective flows of heavy-ion collisions indicate a soft 
EOS [17-21]. In addition, the ratios of the isospin parti- 
cles in heavy-ion collisions, such as 7~/7+, K°/K*, and 
=~ /*, are thought to be sensitive to the stiffness of the 
EOS [22-28]. The production of pions and kaons has been 
experimentally measured in '°7Au+!97 Au collisions. The 
K* production predicted by various transport models favors a 
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soft EOS of isospin-symmetric nuclear matter at high baryon 
densities [29-33]. The m~ /r* ratio predicted by various 
transport models is still model-dependent [34-37]. Based 
on the FOPI data for the n~ /r* ratio [38], some results fa- 
vor a stiff symmetry energy (isospin asymmetric part of the 
EOS) [34, 35], whereas others imply a soft symmetry en- 
ergy [36, 37]. Recently, by analyzing the ratios of charged 
pions in 497Sn +124 Sn, 12Sn +124 Sn, and 19Sn +"? Sn 
collisions [39], a slope of the symmetry energy ranging from 
42 to 117 MeV was suggested [40, 41]. The collective flows 
and ratios of charged pions are still worth studying to find the 
sources of the difference in various transport models and to 
extract information about the EOS from heavy-ion collisions. 


Quantum molecular dynamics (QMD) is a popular trans- 
port model that has been developed into many versions 
and used to successfully describe the properties of nu- 
clear matter, nuclei, mesons, and hyperons [33, 34, 42-55]. 
In high-energy heavy-ion collisions, the relativistic effects 
should be considered in QMD because they become signif- 
icant. The RQMD approach has been proposed for this pur- 
pose [42, 43]. Recently, relativistic mean field theory (RMF) 
was implemented in RQMD (RQMD.RMF) [44-46]. The 
RQMD.RMEF has been successfully applied to investigate the 
collective flows of hadrons [44-46]. In this study, we im- 
plemented RMF theory with isovector-vector and isovector- 
scalar fields in the Lanzhou quantum molecular dynamics 
model (LQMD.RMF). The channels for the generation and 
decay of resonances (A(1232), N*(1440), N*(1535), etc.), 
hyperons, and mesons were included in the previous LQMD 
model [33, 34, 47-50]. Using the LQMD.RMF, we explored 
the relationship between the symmetry energy and the prop- 
erties of the collective flow and pion production. 


The remainder of this paper is organized as follows. In 
Sec. II, we briefly introduce the formulas and approaches 
used in this study. The formulas include RMF theory, the 
dispersion relation, and the production of pions. The results 
and discussion are presented in Sec. III. Finally, a summary 
is presented in Sec. IV. 


72 Il. FORMALISM 


73 A. Relativistic mean field approach 


74 The RMF interaction is achieved by exchanging mesons. 
7s Scalar and vector mesons provide medium-range attrac- 
7s tion and short-range repulsion between nucleons, respec- 
77 tively [56]. The nonlinear self-interaction of the o meson is 
7s introduced to reduce the incompressibility of a reasonable do- 
main [57]. To investigate the properties of symmetry energy, 
so we also consider the isovector-vector p [58] and isovector- 
ai scalar ô mesons [59]. The Lagrangian density is expressed 
s2 as [59, 60] 
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s7 where My = 938 MeV is the nucleon mass in free space, gi 
s With i = 0, w, p, ô is the coupling constant between nucleons, 
s9 M; with i = g, w, p, ô denotes the meson mass, gz and g3 are 
the coupling constants for the nonlinear self-interaction of the 
e1 o meson, and Fy = pwy — O,W,, and Bu = aby — ðb, 
ə2 are the strength tensors of the w and p mesons, respectively. 
93 The equations of motion for the nucleon and meson are ob- 
ə tained from the Euler-Lagrange equations and are written as: 
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98 mo + g20° + 930° = goby = gops (3) 
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mibo = gpI T3 Y = pps; (5) 
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105 Where p and pg are the baryon and scalar densities, respec- 
1% tively, p3 = Pp — Pn is the difference between the proton 
1077 and neutron densities, and ps3 = Psp — Psn is the difference 
108 between the proton and neutron scalar densities. 

In the RMF approximation, the energy density is given by 


m553 = gE PT3 = gspss, 
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where pp denotes the nucleon Fermi momentum, and MX 
Mn — goo F gs563 (— proton, + neutron) denotes the ef- 
fective nucleon mass. The isospin splitting of the effective 
nucleon mass M} — M% still has a large uncertainty at this 
point. Analyses of nucleon-nucleus scattering data based on 
the optical model favor M; — M; > 0 [61, 62]. Calculations 
based on Brueckner theory [63—65] and density-dependent 
relativistic Hartree-Fock [66] also indicate M% — M% > 0. 
However, Mý — M} < 0 is predicted by the transport model 
for heavy-ion collisions [67, 68] and nonlinear RMF mod- 
els [59, 60, 69]. In addition, both M% — Me < 0 and 
M,, — M; > 0 can be predicted by the point-coupling RMF 
[69] and Skyrme and Gogny forces [70-76]. Because the La- 
grangian density in this study is the same as that in Ref. [59] 
and [60] except for the parameter settings, as shown in Fig. 1, 
the negative isospin splitting of the effective nucleon mass 
Mà — Mọ% < 0 is consistent with that in Ref. [59] and [60]. 
In the nonlinear RMF model, the isospin splitting of the ef- 
fective nucleon mass is primarily determined by the coupling 
strength of the 6 meson. The large coupling strength of the 
ô meson results in large isospin splitting of the effective nu- 
cleon mass. When the coupling strength of the ô meson is 
zero (setl), there is no isospin splitting of the effective nu- 
cleon mass. 
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Fig. 1. (Color online) Difference between neutron and proton effec- 
i tive masses as a function of the baryon density. 
137 


138 Using the isospin asymmetry parameter a = (Pn — 
139 Pp)/(Pn + Pp), the symmetry energy can be written as [59] 


1 0° E(p, a) 1 0?[e(p, a)/p] 
140 Esym= 5) la=0 = 2 la=0 
2 Oa 2 Oa 
lpp 1 fs M**p 
141 aE f z , (8) 
6 Eh 2°? 2 EX (1+ fsA(pr, M*)] 
142 Where fi = si = p, 8, and E% = y pp + M*? and M* = 


143 My — goo are the effective nucleon masses of the symmetric 


144 nuclear matter. The integral A(pr, M*) is defined as 
4 p? 
A(pr, M*)= d’ 
(pr, M*) (on I PRT MAA 
_ a PS P 
= M* E% } 


146 (9) 
147 Based on the symmetry energy Esym, the slope L and curva- 
ture Ksym Of the symmetry energy at the saturation density 
po can be obtained as 
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In this study, we set the saturation density to po 
0.16 fm, and the binding energy per particle of the sym- 
metric nuclear matter was set to E/A— My = —16 MeV. For 
symmetric nuclear matter, we set setl, set2, and set3 models 
to be the same as the result of vanishing isospin asymme- 
try. As shown in Table | and Fig. 2, the symmetry energy 
of setl, set2, and set3 was set to be 31.6 MeV at the satu- 
ration density. Setl contained only p mesons; however, set2 
and set3 contained both the p and ô mesons. For setl, when 
the symmetry energy was set to be 31.6 MeV at the saturation 
density, the coupling parameter g, was fixed and the slope 
of symmetry was fixed at L = 85.3 MeV. For set2 and set3, 
when the symmetry energy was fixed at 31.6 MeV, the slope 
of the symmetry energy was obtained by varying the coupling 
parameters gp and gs. Because the effective mass M* of the 
above models for symmetric nuclear matter was the same, the 
symmetry energy with both the p and 6 mesons could not be 
softer than that of set] containing only p mesons. To broaden 
the range of the slope parameters, we set the slope parame- 
ter of set2 and set3 to be 109.3 and 145.0 MeV by varying 
the coupling parameters g, and gs, respectively. A broader 
range of slope parameters would be helpful for understanding 
the relationship between the properties of symmetry energy 
and the observables of heavy-ion collisions. The curvature 
of the symmetry energy Kym, which is a higher-order expan- 
sion coefficient of the symmetry energy compared to the slope 
parameter L, may also affect the properties of the symmetry 
energy and the observables of heavy-ion collisions at densi- 
ties far beyond the saturation density. Ksym of setl, set2, and 
set3 is obtained as -15, 141, and 391 MeV, respectively. 
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182 B. Relativistic quantum molecular dynamics approach 


183 To investigate high-energy heavy-ion collisions, RQMD 
was proposed [42, 43]. Recently, RMF was implemented in 
RQMD [44—46]. In RQMD, for an N-body system, there 
are 4N position coordinates q! and 4N momentum coordi- 
nates pi’ (i = 1,..., N). However, the physical trajectories 
(q; and pi) are 6N for an N-body system. 2N constraints 
189 are required to reduce the number of dimensions from 8N to 
190 physical trajectories 6N [42—46, 77-79]. 


di © O(é = 1,...,2N), 


18. 


A 


18: 


a 


186 
187 


188 


o 


(1) 


191 


150 


set1(L=85.3MeV) 
set2(L=109.3MeV) 
- — — set3(L=145.0MeV) 


100 


p/P, 


Fig. 2. (Color online) Symmetry energy as a function of the baryon 
density. 


where 2N constraints satisfy the physical 6N phase space. 
The sign ~ indicates Dirac’s weak equality. The on-mass 
shell conditions can reduce the phase space from 8N to 7N 
dimensions. 


bi = pi? — M? = (pi — Vi)? — (Mn — Si}? = 0, (12) 


where 2 = 1,..., N. The remaining N constraints are time fix- 
ation constraints. A simple choice of time fixation constraints 
that obey the worldline condition is written as [43, 77, 79, 80] 


ditn =4-(G—gn) =0,( =1,...,N—1), 
2N =G4-qn—T=0, (13) 
where â = (1, 0) denotes a four-dimensional unit-vector [42— 
203 46, 77]. In a two-body center-of-mass system, â is defined as 
2 pi',/\/ pz, with pi’, = p} + pi. We observe that only the con- 
straint 2 = 2N depends on 7. With the above 2N constraints, 
the number of dimensions from 8N will reduce to 6N. These 
2N constraints are conserved over time. 
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2N 
> CRAs = 0. (14) 
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Because only the constraint 7 = 2N depends on 7, A is writ- 
ten as[77] 


a 2N — 1), (15) 


212 with Om = [¢;, d;]. Following previous studies, the Hamil- 
213 tonian of the N-body system was constructed as a linear com- 


214 bination of 2N — 1 constraints [77, 79, 80]: 


2N—1 


H= MSs, (16) 
i=l 
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Table 1. Parameter sets for RMF. The saturation density po is set to 0.16 fm~*. The binding energy of the saturation density is E/A— Mn = 
—16 MeV. The isoscalar-vector w and isovector-vector p masses are fixed at their physical values, mu = 783 MeV and mp = 763 MeV, 
respectively. The remaining meson masses, m, and ms, are set to be 550 and 500 MeV, respectively. 


model] go | gw | g2(fm=")) g3 Ip | 95 | KMeV)| Esym(po)(MeV)} L(oo)(MeV)) M*(po)/My | Ksym (MeV) 
setl | 8.145) 7.570} 31.820 | 28.100) 4.049 - 230 31.6 85.3 0.81 -15 
set2 | 8.145) 7.570} 31.820 | 28.100} 8.673 | 5.347 230 31.6 109.3 0.81 141 
set3 | 8.145) 7.570) 31.820 | 28.100} 11.768| 7.752 230 31.6 145.0 0.81 391 
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Assuming [¢;,¢;] = 0, A, = O for N+1 <i < 2N [77]. 


217 The equations of motion are then obtained as 
N 
dqi 06; 
218 ar T el = Ngy 
N 
dpi 09; 
219 77 = pil = rr (17) 


with the on-mass shell conditions (Eq. (12)) as inputs, the 
equations of motion can be obtained as 
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_ pt M*OM* og, OV; 
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223 Di = =A OF eee Oi: 
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where z;" = p;"/pt° and M* = My — Sj. The scalar 

potential S; and vector potential V;,, in RQMD are defined as 


Si = Jogi + gstidi, 


Visu = Biguri, p ag Bitigpbi,u, (19) 


where t; = 1 for protons, and t; = —1 for neutrons, and B; 
is the baryon number of the ith particle. The meson field can 
be obtained from the RMF: 


mo; + 9207 + 930? = gops., 

miu = Judi, 

m6; = 95(PSp,i — PSn,i) = J5P33,i» 

mb! = Jp(Pp.i — Pn.i) = 9p RY. (20) 


Substituting Eq. (20) into Eq. (19), the scalar potential S; 
and vector potential V;,, of the nucleons can be obtained. 
For other hadrons, such as A resonances, similar to other 
transport models [34, 85], the A optical potential is esti- 
mated using the nucleon potentials and square of the Cleb- 
sch—Gordan coefficient. In the RQMD approach, the scalar 
2a1 density, isovector-scalar density, baryon current, and isovetor 
242 baryon current are expressed as 


M; M; 
243 Psi = 5 =T Pijs Ps3, = X tj pi, 
j#i “Í jai Pi 
p” p” 
244 Jt => Bj pij, R} =X B43. CD 
Era Pj a dD; 
jżi jżi J 


245 Because the difference between the numerical results ob- 
zs tained using the effective mass M* and kinetic momentum 
27 p“ in the density and current and those obtained using a 
2s free mass M; = My 938 MeV and canonical mo- 
249 mentum ph in the density and current is small, a free mass 
M; = My = 938 MeV and canonical momentum ph have 
been used in the above density and current [45]. The interac- 
tion density p;; is given by a Gaussian 
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Pii = CagLg 3/2? ALG 
dijo (P7 +PF) 
y (pitp;)? 
and Ņ;;j is a Lorentz factor that ensures the correct normaliza- 
tion of the Gaussian [81] and is equal to (p? + p?)/(pi + p3) 
in the two-body center-of-mass frame. In this study, we set 

the Gaussian width to Lg = 2.0fm?. 


where dij =q} -| |? is a distance squared [77], 


C. Dispersion relation and production of pions 


The Hamiltonian of the mesons is defined as [48, 82-84] 
Nu 
Hy = ye + w(Di, pi], 


i=l 


(23) 


where V,© is the Coulomb potential, which is expressed as 


a Exe 
C itj 
V.~“ = = 


i 
Tij 


269 (24) 

j=1 
and Nm and Np are the total numbers of mesons and 
25 baryons, including charged resonances, respectively. The 
26 pion potential in the medium, which contains isoscalar and 


isovector contributions, is defined as 


264 


w (Di, pi) = Wisoscalar (Dis pi) + C,720(p/po)”, (25) 
where a denotes the isospin asymmetry parameter, the coef- 
ficient Cr is 36 MeV, the isospin quantity 7 is 1, 0, and -1 for 
na, 7°, and 1+, respectively, and yr determines the isospin 
splitting of the pion potential and is set to two. In this study, 
the isoscalar part of pion potential wisoscalar Was Chosen as 
274 the A-hole model. The pion potential, which contains a pion 
27s branch (smaller value) and A-hole (larger value) branch, is 


276 defined as 
277 


(Dis Pi)Wr—like (Di, Pi) + 
Ji, Pi)WA-—like (Di, Pi): 


Wisoscalar (0 5 pi)= S 
Sa( (26) 
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The probabilities of the pion and A-hole branches satisfy the 
following equation: 


The probability of both the pion and A-hole branches is de- 
fined as [84] 


(27) 


1 
~ 1— OT (w) /dw?’ 
where w denotes wr—like and WA-jike. The eigenvalues of 
Wa—like aNd WA -like are generated from the pion dispersion 
relation 


S (Pi, pi) (28) 


w =p + mÈ +w), (29) 


where II denotes the pion self-energy. Including the short- i 


range A-hole interaction, the pion self-energy is defined as 

ace 
where m, denotes pion mass. The Migdal parameter, g’, was 
set to 0.6. x is defined as 


(30) 


3 
wae exp(—298/02), GD 
A Ww 
where wa = \/m +p; — Mn and ma is the delta mass. 
In this study, the mN A coupling constant fa was 2, and the 
cutoff factor b was TMz. 

Both the Coulomb and pion potentials contribute to the de- 
cay of resonances and the reabsorption of pions. For instance, 
the energy balance of A° in the decay of resonances and the 
reabsorption of pions can be written as 


ye 


Tw) 


(32) 


m3, + = y M8 + (Pr — Be)? + wa (Pr p) + 


where ppr and p, are the momenta of the resonances and pi- 
ons, respectively, and mp is the mass of the resonances. Be- 
cause A° is uncharged, the Coulomb potential exists only for 
charged particles after the decay of A°. 

The pion is generated from direct nucleon-nucleon colli- 
sion and decay of the resonances A(1232) and N*(1440). 
The reaction channels of the resonances and pions, which 
are the same as those in the LQMD model, are as fol- 
lows [33, 48, 88, 89]: 


NNGNA, NNGNN*, NNGAA, ASN, 
N* & Nr, NN > NNn(s — state). (33) 


For the production of the A(1232) and N* (1440) resonances 
in nucleon-nucleon scattering, the parameterized cross- 
sections calculated using the one-boson exchange model were 
employed [91]. 

The decay width of A(1232) and N*(1440), which orig- 
inates from the p-wave resonances, is momentum-dependent 
and is expressed as [91] 


a.\pl° 
(1 + a2|p1?)(as + pI?) 


T(|pl) = To, 


(34) 
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where |p] is the momentum of the created pion (in GeV/c). 
The parameters a1, a2, and a3 were taken as 22.48 (17.22) 
c/GeV, 39.69 (39.69)c?/GeV?, and 0.04(0.09) GeV?/c?, 
respectively, for A(N*). The bare decay width of A(N*) 
was given by T'o = 0.12(0.2)GeV. With the momentum- 
dependent decay width, the cross-section of pion-nucleon 
scattering has the Breit-Wigner form: 
0.251? (P) 


Onn (Vs) = Omax( Trp + (Vs — mo) 


where p and pm are the three momenta of the pions at ener- 
gies of \/s and mo, respectively. The maximum cross-section 
Omax Of A and N* resonances was obtained by fitting the 
total cross-sections of the experimental data in pion-nucleon 
scattering using the Breit-Wigner formula [92]. For instance, 
the maximum cross-section Gmax of A resonance was 200, 
133.33, and 66.7 mb for rtp + Att (mn > A), 
np + At (mn > A®) and a-~p— A} (ntn > AT), 
respectively [89]. 

Note that the threshold effect was neglected in this 
study. The threshold effect mainly refers to the A pro- 
duction threshold energy and incident energy of two col- 
liding nucleons modified by the medium. The inci- 
dent and threshold energies in the medium are defined 


as A/S = yE EE? + Eš +59)? (Sy + 2)? and 
Jah = Vimy + 9 + mä + D2)? — (U3 + Ea), respec- 
tively [85-87]. Because Š; = 0 and p} ~ 0 for static nuclear 
matter, the difference between the incident and threshold en- 
ergies is \/Sin — Sin = Ef + X? + ER + X9 -— m3 — 
£9 — mi — X} [86]. The difference between the incident and 
threshold energies, which is isospin dependent, may result in 


an enhancement in the nn — pA- channel and suppression 
of the pp — nA** channel. 


Pm 
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Il. RESULTS AND DISCUSSIONS 


The directed and elliptic flows were derived from the 
Fourier expansion of the azimuthal distribution: 


N pr) = Noll +2Vi(y,pr)cos(¢) 


do 
+2V2(y, pr )cos(2¢)], (36) 
where the azimuthal angle of the emitted particle ø was mea- 
sured from the reaction plane. pr = ,/p? + ite denotes the 


transverse momentum, and the directed flow V; and elliptic 
flow V2 are expressed as follows: 


Sao 2s 
PT 
2, 
Py 
The directed flow provides information on the azimuthal 
anisotropy of the transverse emission. The elliptic flow de- 
scribes the competition between the in-plane (V> > 0) and 

out-of-plane (V2 < 0) emissions. 


Vı =< cos(¢) 


V2 =< cos(2¢) >=< Su (37) 
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Firstly, the collective flows of LQMD.RMF in the 
197 Au +197 Au collisions have been investigated at an 
incident energy of 2.92 GeV/nucleon (the corresponding 
nucleon-nucleon center-of-mass energy is /Svy = 3 GeV 
). The collective flows of LQMD with Skyrme interaction and 
without the momentum-dependent interaction have also been 
investigated. The Skyrme interaction of symmetric nuclear 
matter is taken to be the same as that of SLy6, with an incom- 
pressibility of K = 230 MeV at po = 0.16fm~°[93, 94]. 
The symmetry energy of Skyrme interaction is defined as 

2 

Esym = 3 Ea (3 2)2/8 a Coym (p/ po)”. When the Csym 
and ys; are set to be 38.6 MeV and 1.049, respectively, the 
symmetry energy and the slope parameter of symmetry en- 
ergy are 31.6 MeV and 85.3 MeV, respectively. The symme- 
try energy and the slope parameter of symmetry energy of the 
Skyrme interaction are as same as those of setl. As shown in 
Fig.3, we have compared the collective flows of LQMD.RMF 
and LQMD with Skyrme interaction with recent experimen- 
tal data from STAR Collaboration[90]. The collective follows 
of LQMD.RMF and LQMD with Skyrme interaction can de- 
scribe the STAR data at an impact parameter b=4 fm and b=7 
fm, respectively. The directed flows of LQMD.RMF are al- 
most consistent with the STAR data across the entire rapidity. 
However, The directed flows of LQMD with Skyrme interac- 
tion are weaker than the STAR data across the entire rapid- 
ity. This result may be due to the fact that the value of di- 
rected transverse momentum with the Lorentz effect is larger 
than that without the Lorentz effect[43]. The elliptic flows of 
both LQMD.RMF and LQMD with Skyrme interaction are 
consistent with the STAR data at rapidities smaller than 0.5. 
However, the elliptic flows of both LQMD.RMEF and LQMD 
with Skyrme interaction are weaker than the STAR data at 
large rapidity. At an incident energy of 2.92 GeV/nucleon, the 
LQMD.RMEF can better describe the experimental data com- 
pared to the LQMD with Skyrme interaction and without the 
momentum-dependent interaction. Based on the above analy- 
ses, the RMF models have been implemented into the LQMD 
model successfully. 


With this LQMD.RMF model, the !°8Sn +112 Sn and 
1325n +174 Sn collisions in this study were investigated at 
an incident energy of 270 A MeV and an impact parameter 
b=3 fm. At an incident energy of 270 A MeV, the nuclear 
matter of the collision center can be compressed to densi- 
ties approaching 2. In this dense region, collective flows, 
which reflect nucleon-nucleon interactions, can be used to ex- 
tract the high-density behavior of the EOS [1, 44, 79, 88, 95]. 
The collective flows of free protons in the !9°Sn +1! Sn and 
132Sn + !24Sn collisions are shown in Fig. 4 and Fig. 5, re- 
spectively. It is reasonable that the maximum value of the 
directed flow V; was significantly larger than that of the el- 
liptic flow V2. In the same reaction system, the difference 
in the directed flows with various slopes of symmetry energy 
(setl, set2, and set3) was small. The difference in the ellip- 
tic flows with various slopes of symmetry energy was also 
small. To determine the relationship between the slope of the 
symmetry energy and the collective flow, we must process the 
collective flow data. The difference between neutron and 
proton directed flows emitted from heavy-ion collisions can 
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Fig. 3. (Color online) Rapidity distribution of the collective flows 
of free protons in the 1°’ Au +'°" Au reaction at an incident energy 
of 2.92 GeV/nucleon. The open circles correspond to the STAR 
data[90]. 


be used to extract the density dependence of the symmetry 
energy [88, 95]. The difference between the neutron and pro- 
ton directed flows is defined as Vin — Vip, as shown in Fig. 6. 
The trend and shape of the difference between the neutron and 
proton directed flows were similar to those of nonrelativistic 
LQMD [88]. For a given reaction system (nearly symmetric 
108Sn +112 Sn system or neutron-rich !8*Sn +124 Sn system), 
the absolute value of the difference between the neutron and 
proton directed flows with a soft symmetry energy was higher 
than that with a stiff symmetry energy. Interestingly, the stiff- 
ness of the symmetry energy can be reflected through the dif- 
ference between the neutron and proton directed flows. The 
relationship between the curvature of the symmetry energy 
K sym and the collective flows was then investigated. The 
difference between K gym of setl and Ksym of set3 was 406 
MeV, which was significantly larger than the 59.7 MeV differ- 
ence between L of set! and L of set3. Although the curvature 
of the symmetry energy Ks, also affected the difference be- 
tween the neutron and proton directed flows, because the nu- 
clear matter of the collision center could only approach 2p at 
an incident energy of 270 A MeV, the effects of Ksym were 
not significant compared to the effects of the slope parameter 
L. 

In addition to collective flows, the production of isospin 
exotic particles such as hyperons, kaons, and pions can also 
be used to extract the symmetry energy [22-28]. Because the 
thresholds of hyperons and kaons were not reached at the in- 
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Fig. 4. (Color online) Rapidity distribution of the collective flows of 
free protons in the ‘°°Sn +"? Sn reaction at an incident energy of 
270 MeV/nucleon. 
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Fig. 5. (Color online) Rapidity distribution of the collective flows of 
free protons in the ®*Sn +14 Sn reaction at an incident energy of 
270 MeV/nucleon. 


458 
459 
460 
46 


462 
463 
464 
465 


466 


0.04 
18Sh+"Sn 
F 270MeV/nucleon 
0.02 set1(L=85.3MeV) 
Su —-—--set2(L=109.3MeV) 
ee — — — set3(L=145.0MeV) 
0.00 
‘ce -0.02 
> 
0.00 
-0.04 1325+ Sn 
270MeV/nucleon 
-1.0 -0.5 0.0 0.5 1.0 
yly 


proj 
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an incident energy of 270 MeV/nucleon. 
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Fig. 7. (Color online) Ratio between the 7~ and 7* yields as a 
function of the incident energy in the '°*Sn +14 Sn reaction. 


cident energies in this study, the isospin exotic particles were 
pions. First, we calculated the ratio between the 7~ and at 
yields of the neutron-rich 1°?Sn +124 Sn system as a func- 
tion of the collision energy at the impact parameter b=3 fm 
and Oem < 90°. Because setl had the softest symmetry en- 
ergy, it had the highest neutron density. Consequently, there 
were more neutron-neutron scatterings in setl, resulting in 
the production of more A~ and 7~. As shown in Fig. 7, the 
a /m* ratio of setl was the highest, and the m~ /7* ratio 
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Fig. 8. (Color online) Transverse momentum spectra of pions as 
functions of transverse momentum at an incident energy of 270 
MeV/nucleon. The left two panels [(a) and (b)] are the results of 
the 1°°Sn +11? Sn reaction, and the right two panels [(c) and (d)] 
are the results of the 1??Sn +174 Sn reaction. The full circles corre- 
spond to the S7RIT data [40]. 


of set2 was higher than that of set3. Specifically, at a colli- 
sion energy of 270 MeV/nucleon, the m~ /7* ratio without 
(with) the r potential changed from 2.71 (2.54) to 2.23 (2.06) 
when the slope parameter was varied from L =85.3 to 145.0 
MeV, that is, from setl to set3. In other words, the 17 / at 
ratio as a function of collision energy depends on the stiff- 
ness of the symmetry energy. This result was consistent with 
the predictions of most transport models [28, 39, 86]. When 
the r potential was considered, the interaction between 7 and 
the nucleon became attractive, resulting in a decrease in both 
n~ and 7+ via the absorbed channels nN —> A(1232) and 
A(1232)N — NN. However, with the 7 potential, because 
there were more neutrons in the neutron-rich 82Sn +174 Sn 
system, 7~ was more easily absorbed than 7*. Consequently, 
the m~ /r* ratio without the 7 potential was higher than that 
with the 7 potential in the same RMF model. Moreover, it 
is worth mentioning that the threshold effect neglected in this 
study may cause the m~ /r* ratio to be reversed. In other 
words, with the threshold effect [85-87], m~ /7* of set3 may 
be the highest, and the 7~/z* ratio of set2 may be higher 
than that of set1. 

Next, the properties of 7 were predicted as a function of 
the transverse momentum. As shown in Fig. 8, the left and 
right panels are the transverse momentum spectra of pions 
for the nearly symmetric 1°°Sn +!!? Sn and neutron rich 
132Sn +14 Sn reactions at ĝem < 90°, respectively. In col- 
lisions between isotopes, 7* is mainly generated from colli- 
sions between protons and m~ is mainly generated from col- 
lisions between neutrons. Theoretically, a stiffer symmetry 
energy would have a stronger repulsive force to push out neu- 
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Fig. 9. (Color online) Single spectral ratios of pions as functions 
of transverse momentum for the !°°Sn +*!? Sn and '8?Sn +14 Sn 
reactions at an incident energy of 270 MeV/nucleon. The full circles 
correspond to the S7RIT data [40]. 


trons and a stronger attractive force to squeeze protons, re- 
sulting in a decrease in the m~ yield and an increase in the 
mt yield, respectively. As shown in panels (b) and (d) of 
Fig. 8, the stiffer symmetry energy indeed led to larger trans- 
verse momentum spectra for 7+. However, the relationship 
between the transverse momentum spectra of m~ and the stiff- 
ness of the symmetry energy could not be directly explained. 
Compared with the stiffness of the symmetry energy, the 7 
potential had a more significant impact on the transverse mo- 
mentum spectrum of 7. For the neutron-rich !8*Sn +174 Sn 
system, the transverse momentum spectra of both 7+ and m7 
without the 7 potential were lower than those of the StRIT 
data at pr < 200 MeV. When the m potential was consid- 
ered, the transverse momentum spectra of both a? and m7 
increased at pr < 200 MeV but decreased at pr = 200 MeV. 
Consequently, the transverse momentum spectra of ++ with 
the m potential were almost consistent with the S7RIT data 
[40]; however, the transverse momentum spectra of 7~ were 
still lower than the S7RIT data for the entire pp domain. The 
lower transverse momentum spectra of 7~ may be due to the 
absence of a threshold effect. The threshold effect, which was 
not considered in this study, may enhance the production of 
n~ [85-87]. 

Because a stiffer symmetry energy would have a stronger 
repulsive force to push out neutrons and a stronger attractive 
force to squeeze protons, resulting in a decrease in the m~ 
yield and an increase in the 7* yield, respectively, the single 
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Fig. 10. (Color online) Double ratio of pions as a function of trans- 
verse momentum at an incident energy of 270 MeV/nucleon. The 
full circles correspond to the S7RIT data [40]. 


ratio SR(n—/n*) = [dM (x—)/dpr]/[dM(x*)/dpr] may 
depend on the stiffness of the symmetry energy and the reac- 
tion system. As shown in Fig. 9, for both the nearly symmet- 
ric system and neutron-rich system, a softer symmetry energy 
led to a larger single ratio. In addition, for the same stiffness 
of the symmetry energy, because the neutron-neutron scatter- 
ing of the neutron-rich !32Sn +174 Sn system was greater than 
that of the nearly symmetric }°°Sn +11? Sn system, the sin- 
gle ratio of the neutron-rich system was higher than that of 
the nearly symmetric system. It is worth mentioning that the 
single ratio of °°Sn +11? Sn was almost consistent with the 
experimental data. However, the single ratio of 13?Sn +1°4 Sn 
was lower than that of the experimental data for the entire pr 
domain. This was because the transverse momentum spectra 
nT of 18?Sn +124 Sn were lower than the experimental data. 


The double ratio between the neutron rich system 
and the nearly symmetric system DR(zx~/77*) 
SR(m7 /a*)1324124/SR(a7 /7*) 1084112; which can 
cancel out most of the systematic errors caused by Coulomb 
and isoscalar interactions, was suggested to extract the 
properties of the symmetry energy [40]. However, because 
a lower symmetry energy will lead to a larger single ratio for 
both the nearly symmetric system and neutron-rich system, 
as shown in Fig. 10, it is still difficult to understand the 
dependence of the double ratio on the stiffness of symmetry 
energy. In addition, the double ratio without the m potential 
decreased slightly as a function of the transverse momentum, 
whereas it increased slightly as the transverse momentum 
increased when the 7 potential was considered. However, the 
increasing trend was still considerably weaker than that of 
the experimental results. The lower double ratio originated 
from the lower single ratio of the neutron-rich 1°?Sn +174 Sn 
system compared with the experimental data. The threshold 
effect may enhance the production of m~ [85-87] and the 
single ratio of the neutron-rich system, resulting in a larger 
double ratio. 
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IV. CONCLUSION 


An RMF with various symmetry energies, namely setl, 
set2, and set3, was implemented in LQMD. The collective 
flows of the nearly symmetric }°°Sn +1!? Sn and neutron- 
rich '?Sn +14 Sn systems were successfully generated from 
the LQMD.RMF. It has been observed that the directed flow 
Vı was an order of magnitude larger than the elliptic flow 
Və. However, the difference between the directed flows V; 
with various slopes of symmetry energy was small. The dif- 
ference in the directed flows V2 with various slopes of sym- 
metry energy was also small. To explore the relationship be- 
tween the collective flow and the stiffness of the symmetry 
energy, we defined the difference between neutron and proton 
directed flows Vin — Vip. For a given nearly symmetric sys- 
tem or neutron-rich system, the absolute value of Vin — Vip 
increased with decreasing slope of the symmetry energy. 


We also investigated the relationship between isospin ex- 
otic particles and the stiffness of the symmetry energy. The 
ratio between m~ yield and m* yield of the neutron-rich 
132Sn +124 Sn system as a function of the collision energy 
increased with a decrease in the slope parameter of the sym- 
metry energy. At an incident energy of 270 MeV/nucleon, the 
properties of 7 were predicted as a function of the transverse 
momentum. For the nearly symmetric !°°Sn +11? Sn system, 
the single ratio of the nearly symmetric system was consistent 
with the experimental data. However, for the neutron-rich 
1326n +124 Sn system, the single ratio was lower than the 
experimental data, resulting in a double ratio lower than the 
experimental data. The 7 potential did not explain the lower 
transverse momentum spectra of 7~ in the neutron-rich sys- 
tem. Considering the m potential, the double ratio increased 
slightly with increasing transverse momentum. However, the 
increasing trend was still considerably weaker than that ob- 
served in the experimental results. The single ratio of the 
neutron-rich system and the double ratio may be lower than 
the experimental data because of the lack of a threshold ef- 
fect. The threshold effect, which can enhance the production 
of m~, could be a candidate for enhancing the single ratio of 
the neutron-rich system to a double ratio. Moreover, because 
a softer symmetry energy led to a larger single ratio for both 
nearly symmetric and neutron-rich systems, the dependence 
of the double ratio on the stiffness of the symmetry energy 
was not significant. The sensitivity of the double ratio to the 
stiffness of the symmetry energy may also be due to the lack 
of a threshold effect. When the threshold effect is considered, 
the production of 7~ in a neutron-rich system may be more 
significant than that in a nearly symmetric system. In the fu- 
ture, when collective flows occur, and the single ratio of the 
neutron-rich system and the double ratio of the LQMD.RMF 
are consistent with the experimental data, Vin — Vi, and the 
single ratio of the neutron-rich system m~ /7*, which are sen- 
sitive to the stiffness of the symmetry energy, may be used to 
extract the slope parameter of the symmetry energy. 
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Appendix A: DETAILS of EQUATION OF MOTION 


For numerical calculations, the equation of motion (Eq. 
(18)) must be written in the computed form. With Eq. (19) 
and M¥ = Mi — S; = My — S; as the inputs, the equation 
of motion (Eq. (18)) can be expanded as 
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where z! = p} /p® and z;" = p;"/p*°. Based on Eq. (20), 


poe and Ta are obtained as follows: 
00; 06; 
Oi = 5 Go z? (3 = s (A8) 
Opsi mZ +2g20; + 39307 ps3; = MF 


In the two-body center-of-mass frame, p;; equals p;;. The 
. 2 2 —_@ CaP) 
squared distance q7 ; jis reduced to q7 ; j= ae 
In the actual calculation, we replace p? with \/p? + M? to 
save calculation time [79]. Thus, the partial derivative of den- 
sity versus momentum and space can be written as 
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where (3;; is defined as (p; + p)/(p? + p$). In addition, 35 
can be written as -Hf The partial derivative of the energy 
0 
component of Dag is zero and that of the momentum compo- 
. =) 
nent of al is written as sit: Using the above equations, the 
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momentum and space of neutrons and protons are known. 
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